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Quark back reaction to deconfinement transition via gluon propagators 
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Contribution of the quark back reaction to the deconfinement phase transition is studied in the 
thermodynamical potential of the gluonic sector which consists of the gluon and ghost propagators 
calculated in the lattice QCD simulations. Starting from QCD thermodynamic potential, we define 
the gluonic potential in the leading-order of the 2PI formalism, which can describe the deconfinement 
phase transition. Then the gluonic potential can be written by using the microscopic characters; 
the gluon and ghost propagators in Landau gauge fixing. To include effects of the quark back 
reaction, we calculate the gluon propagators in lattice QCD simulations with two-flavored dynamical 
quarks. Fitting the lattice data by the Gribov-Stingl form and investigating the phase transition 
of the gluonic potential, we find that enhancement of the quark back reaction reduces the critical 
temperature of the deconfinement phase transition. 
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I. INTRODUCTION 

Understanding of a nature of the deconfinement phase 
transition in Quantum Chromodynamics (QCD) is a long 
standing problem in the nuclear physics [1]. Experi- 
mental seeking of the phase transition and quark-gluon 
plasma state has been done by the relativistic heavy- 
ion collision in RHIC and LHC [2]. To reveal nontriv- 
ial properties of the phase transition, it is important to 
know thermodynamics of quarks and gluons where non- 
perturbative aspects of QCD play important roles. 

A lattice QCD simulation is a powerful tool to deal 
with the non-perturbative effects of QCD and obtain 
computational data on the basis of the first principle 
calculations. On the other hand, an effective theory is 
a useful approach to extract physical meanings and pic- 
tures from the numerical data. Establishment of a re- 
liable effective approach leads us to deep understanding 
of the deconfinement phase transition at finite tempera- 
ture. So far, several approaches have been proposed. The 
Nambu-Jona-Lasinio (NJL) model is the most effective 
one to describe thermodynamics of the quark sector in 
terms of breaking and restoration of the chiral symmetry 
at finite temperature [3]. However the NJL model is not 
approachable to describe thermodynamics of the gluonic 
sector which is probably essential to the deconfinement 
phase transition. Recently, an extended NJL model with 
the Polyakov loop, so called PNJL, has been proposed 
[4, 5]. Although PNJL model can describe approximate 
deconfinement transition, it is not desirable model. This 
is set only by macroscopic properties of the pure gauge 
theory (e.g. pressure in quenched lattice QCD), since 
all components of the gluonic potential consists of the 
Polyakov loop. Therefore, there is no definite procedure 
to introduce microscopic contributions of quark excita- 
tions which couple to the gluonic sector, i.e. quark back 
reaction. 

Microscopic description of the effective approach prob- 
ably inspires us with intuitive understanding of the de- 
confinement phase transition. Such an effective approach 



should consist of the gluon and ghost propagators which 
can represent fundamental degrees of freedom of the 
gluonic sector. The gluon and ghost propagators also 
play an important role in so-called Kugo-Ojima/Gribov- 
Zwanziger confinement scenario [6, 7]. So far, the gluon 
and ghost propagators in the Landau gauge fixing at 
zero and finite temperature have been studied in the lat- 
tice simulations, the Dyson-Schwingcr (DS) equation ap- 
proach, and the functional renormalization group (FRG) 
approach [8-10]. Recently, FRG approach has found 
that the deconfinement transition can be described by 
the Polyakov loop via the Landau-gauged gluon and 
ghost propagators [11, 12], and the quark back reaction 
contributes to critical temperature of the phase transi- 
tion [13]. The gluon and ghost propagators in FRG ap- 
proach are well described in the pure-gauge limit. How- 
ever, contributions arising from the quark back reaction 
are expressed by the quark-meson model which starts 
from the bare mesonic potential at ultra-violet scale and 
represents the quantum fluctuation via the meson propa- 
gations. Therefore it is uncertain to describe appropriate 
quark excitations based on the first principle of QCD. 

To establish sound effective approach from the micro- 
scopic degrees of freedom, it has been proposed that the 
gluonic potential is formulated by starting from the pure 
gauge theory based on the leading contributions of the 
2PI formalism [14]. This approach enables us to deal with 
the propagators obtained in the lattice simulations which 
include all information of the non-perturbative aspects 
of QCD. They have found that several physical quan- 
tities (e.g. the critical temperature, equation of states, 
etc.) can be reproduced via the approach with input of 
the propagators determined in pure-gauge lattice simu- 
lations. 

In this paper, we develop the previous approach to join- 
ing the gluon propagators determined in the lattice QCD 
simulations with two-flavored dynamical (sea) quarks at 
finite temperature and investigate the effect of the quark 
back reaction to the deconfinement phase transition. At 
finite temperature, the four-dimensional transverse gluon 
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propagator can be decomposed into the tree-dimensional 
longitudinal and transverse components. We investigate 
temperature dependence of the propagators in lattice 
QCD simulations. By using these propagators with the 
Gribov-Stingl form, we calculate the critical temperature 
via the effective potential and investigate the effect of the 
quark back reaction to the phase transition of the gluonic 
sector. 

This paper is organized as follows. In the next section, 
we introduce the formalism of the effective approach and 
calculate the thermodynamical potential of the gluonic 
sector. The lattice QCD formalism and simulations are 
shown in Sec. Ill and the critical temperature calculated 
by the gluonic potential is shown in Sec. IV. The paper 
is summarized in Sec. V 



with 4d transverse (-D(t)) an d longitudinal (-D(l)) com- 
ponents. Since the Lorentz symmetry is explicitly broken 
at finite temperature, the 4d transverse gluon propagator 
is decomposed into the three-dimensional (3d) transverse 
and longitudinal components: 



T^D {T) (p 2 ) = P^D T + P£ V D L . 



(5) 



where we impose the Landau-gauge fixing. Then the 
gluon propagators can be defined by the gluon fields as 
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II. GLUON AND GHOST EFFECTIVE 
POTENTIAL 

First of all, let us briefly explain formulation of the 
effective potential of the gluonic sector. In general, the 
effective approach begins to distinguish the QCD ther- 
modynamical potential to matter and gluonic sectors, 



QCD 



u. 



(i) 



One of the descriptions of the matter part is, for instance, 
NJL model. Mixture between f2 ma t and U is induced by a 
classical field <f> which appears in the covariant derivative 
of n mat . We can define it as the temporal components 
of the gluon field: 4> = 03^3 + </>sAs where A3 and As 
are diagonal components of Gell-Mann matrices. At zero 
chemical potential, only A3 path becomes relevant, and 
we can rewrite it only with one angular parameter 6: 
4> = 2ir x diag(6>, —9, 0). Then the Polyakov loop can be 
written as 



2cos(2tt0)]. 



(2) 



In this paper, we focus only on the thermodynamical 
potential of the gluonic sector and construct it by the 
gluon and ghost propagators in the Landau gauge fix- 
ing. The effective potential of the gluonic sector can be 
evaluated by the leading-order contribution of the 2PI 
formalism in the pure-gauge theory [14], 
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tr In D^ 1 , 



(3) 



where "tr" acts on all indices. Quantum fluctuation of <fi 
and spatial gluon fields have already been integrated out 
and contributes via the gluon Da and ghost D c prop- 
agators. Da can be expressed by the four-dimensional 
(4d) transverse (T^ u ) and longitudinal {L^ v ) projection 
tensor, 



Da(p 2 



T^D^t) 



(P 2 ) 



)(/), 



(4) 



where N g = N 2 — 1. In the following, we impose the zero 



Matsubara-frequency on the propagators, i.e. 
The effective potential (3) becomes 



Pi 



0. 
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(8) 



where traces of the Lorentz indices are already done. 

Since analytic forms of the propagators help us to eval- 
uate the effective potential in practice, we assume the 
ghost propagator as 



D C (P 2 ) 



P 



d- 2 



(P 2 ) 2 



(9) 



and the gluon propagator as the "Gribov-Stingl" form: 

c 2 d 2 (p 2 + d- 2 ) 

D (P ) = (p2 + r 2)2 ■ ( 10 ) 

where subscripts L and T for the gluon propagator and 
parameters are omitted. The latter form can be derived 
in a low-energy effective reduction of the Yang-Mills the- 
ory [15]. For the ghost propagator, we utilize results 
of Dc obtained in Ref. [14] by fitting the data in pure- 
gauge lattice simulations [10], since the ghost propagator 
is not sensitive to the dynamical quarks due to lack of 
the quark-ghost vertex [16]. 

By substituting Eqs. (9) and (10), the effective poten- 
tial (8) becomes 

"=(i-2)v pc ,-i(v^(d Z 2 )-2V^(r!)) 

-\(v A T \d- 2 )-2V A T \r 2 T )) +Vc(d c 2 ), (11) 

(L) (T) 

where Vpcri^V^ , V A and Vc are the pcrturbative one- 
loop [17, 18], 3d longitudinal gluon, 3d transverse gluon 
and ghost effective potentials, respectively. The gluon 
and ghost effective potentials can be expressed as 



V(m 2 ;8) 



dp 
2^2 



-Ttr, / ^Wln(l - L 8 e- fj \ /p2+m2 ) , (12) 
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where tr c is the adjoint color trace and L% is the 
Polyakov-loop operator in the adjoint representation 
(L 8 ) a b = 2ti c (t a e^tbe 1 ^). Equation (11) means that, if 
there are the parameters (dj,, rj,, dj>, ...) determined in 
the lattice QCD simulations, one can calculate the glu- 
onic potential from microscopic points of view: The gluon 
and ghost propagators which basically manifest all non- 
perturbative aspects of QCD including the quark back 
reaction. 

It is worth to mention that in the derivation of PNJL 
model, the spatial gluons are expressed by the Polyakov 
loop via the strong coupling expansion and the quantum 
fluctuation of the temporal gluons is suppressed due to 
an application of the classical Haar measure [4]. This 
means that there is no microscopic component (like a 
propagator) in the gluonic sector of PNJL model, and it 
is impossible to deal with the quark back reaction based 
on a definite manner. 



III. GLUON PROPAGATORS ON LATTICE 

In order to investigate effects of the quark back re- 
action to the phase transition of the gluonic sector, it 
is necessary to calculate the gluon propagators in the 
lattice QCD simulations with dynamical quarks. We 
utilize gauge configurations generated by WHOT-QCD 
Collaboration on iV 3 x N t — 16 3 x 4 lattice with a 
rcnormalization-group improved gauge action S g and a 
clover improved Wilson quark action with two flavors S q : 

s 9 = -0E(<* E W £H*) + c i E Kv 2 w) > ( 13 ) 

X fl<V fi^v 

s « = E E^W= ( 14 ) 

/=1,2 x,y 

where c\ = —0.331, c = 1 — 8ci and WJ}* m is a n x m 
shaped Wilson loop. The quark matrix is 

D xy = S xy - if y^[(l - Jp)Ux,n6 x +hy 

(15) 

where K is the hopping parameter, Fa V is the lattice 
field strength and csw is the clover coefficient deter- 
mined in the one- loop perturbation theory: csw — (1 — 
0.8412/3 -1 ) -3 ' 4 . The simulations have been performed 
along the line of constant physics at mpg/my = 0.65 
(ratio of the pseudo-scalar and the vector meson masses 
at zero temperature). Details of the lines of constant 
physics and the phase diagram with the same lattice ac- 
tions are given in Rcfs. [19, 20]. For the investigation 
of the phase transition, we calculate the propagators in 
the temperature range of 0.8 ^ T/T pc £ 1.4, where T pc 
is the pseudo-critical temperature for the transition from 
the hadronic phase to the quark-gluon plasma phase in 
N f = 2 QCD. 



The purpose of the lattice simulations is to evaluate 
the gluon propagators (6) and (7). The gauge field on 
the lattice can be defined by the link valuable as 

where a is a lattice spacing and we assume the isotropic 
lattice. Imposing the Landau gauge fixing on the lat- 
tice configurations, we define the zero-frequency (j?4 = 0) 
gluon propagators in the momentum space, 

9 4=1 X 

D L (p) = ^^tr[A 4 ( 2 ;)A 4 (0)] e " 2 ^^/ JV % (18) 

where kj = 0, 1, • • • , N s — 1 and the physical momentum 
is 

-^X> 2 (^> (19) 

Notice that the transeverse condition PiAi = is imposed 
on D T (p). We measure the propagators with 5000 tra- 
jectories and estimate statistical errors by the jackknife 
method with bin size of 100 trajectories. 

Figure 1 shows temperature dependence of Dt (top) 
and Dl (bottom) as a function of momentum, where a 
physical scale is normalized by T pc ~ 180 MeV. We find 
that there is significant temperature dependence at infra- 
red region (i.e. small p). On the other hand, the prop- 
agators converge at ultra-violet region where the physics 
is insensitive to temperature and probably dominated by 
perturbative descriptions. When temperature increases, 
magnitude of both propagators increases at T < T pc and 
shows enhancement at T ~ T pc for Dl and T ~ 1.2T pc 
for Dt and decreases subsequently. In the case of pure 
gauge theory, the gluon propagators are almost insensi- 
tive to temperature below T c , whereas these monotoni- 
cally decrease above T c [10]. Therefore the strong sensi- 
tivity to temperature below T pc in our case implies that 
there are a large number of configurations which can cou- 
ple to the gluon fields even at T < T pc where the gluon 
confinement is slightly broken by the quark back reac- 
tion. Since the number of such configurations would be 
proportional to magnitude of the gluon propagators, the 
quark back reaction enhances with temperature increases 
at T < T pc . 

Furthermore, the enhancement of the propagators in 
momentum space means strong correlations between glu- 
ons in real space. In general, the longitudinal propagator 
is dominated by the electric gluon sector which is sensi- 
tive to the deconfinement phase transition, e.g. Polyakov 
loop which consists of a sum of the longitudinal gluons. 
It is therefore natural that the correlation of longitudinal 
gluons becomes strong at T ~ T pc . On the other hand, 
the transverse propagator is dominated by the magnetic 
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FIG. 1: Results of the transverse (top) and longitudinal (bot- 
tom) gluon propagators in momentum space for several tem- 
perature. A physical scale is normalized by T pc . 



gluon sector which is probably insensitive to the phase 
transition. For instance, the confinement of spatial Wil- 
son loop, characterized by the spatial string tension, does 
not disappear even above T pc [19]. However the enhance- 
ment of Dt at T ~ 1.2T pc implies that the magnetic sec- 
tor may also play an important role in the deconfinement 
phase transition. 



IV. QUARK BACK REACTION TO EFFECTIVE 
POTENTIAL 

In this section, we calculate the critical temperature 
from the thermodynamic potential of the gluonic sector 
(11) by using the gluon and ghost propagators on the 
lattice, and investigate effects of the quark back reaction 
to the phase transition. To input the lattice results to 
the thermodynamic potential, we fit the gluon propaga- 
tors on the lattice shown in Fig. 1 by the "Gribov-Stingl" 
form (10). Since, in general, a finite size effect becomes 



serious at small p [10], we cut the fit range at some mo- 
mentum, i.e. taking the fit range at p > (p/T pc ) m i n . To 
consider the fit range dependence, we take two values of 
(p/Ipc)mm- We also use only the data at p = (p, 0,0) 
to reduce artifacts of rotational symmetry breaking on 
the lattice. The fit is performed by minimizing x 2 an d 
results are summarized in Tab. I. 

Results show that at T/T pc = 0.83 and 0.86 the fit 
provides good x 2 s, whereas it does not work well when 
temperature increases. The fit range dependence are less 
than 20 % at low temperature, whereas it becomes huge 
at T/Tp C = 1.00. Therefore, we restrict applicable tem- 
perature of the "Gribov-Stingl" form and our effective 
approach within T/T pc < 0.94. 

Then the critical temperature of the effective potential 
T' c can be obtained by the following steps: 

1. Determine the parameters t and at each 
T/T pc . 

2. Minimize the effective potential (11) with respect to 
the classical field </>. 

3. Determine T' c at the point that the Polyakov-loop 
becomes non-zero value. 

Notice that T' c does not correspond to T pc on the lat- 
tice because it is obtained by the purely gluonic sector 
(11). Qualitative shape of the effective potential (11) is 
the same with that of the previous study using the gluon 
propagators in pure gauge theory [14]. Therefore it ex- 
hibits the first order phase transition. Contribution of 
the quark back reaction appears as a shift of the charac- 
teristic scale T' c in the effective approach. 

Table II summarizes results of T' c in GeV unit for sev- 
eral T/T pc . We estimate T' c as an average of two fit ranges 
shown in Tab. I and systematic errors due to the fit range 
dependence as difference from the averaged value, which 
indicate the second parenthesis in the Table. We can see 
that T' c is sensitive to T /T pc and decreases with temper- 
ature increasing at T < T pc . Since, as mentioned above, 
the number of lattice configurations which couple to the 
gluon fields increases with temperature increases at be- 
low T pc , the decrease of T' c is caused by the enhancement 
of the quark back reaction via the gluon propagators. 



V. SUMMARY 

We have investigated effects of the quark back reaction 
to the gluonic thermodynamical potential which consists 
of the gluon and ghost propagators. We have calculated 
the gluon propagators of the 3-dimensional longitudinal 
and traverse sectors in lattice QCD simulations with two- 
flavor dynamical quarks at finite temperature. The sim- 
ulations have been performed on 16 3 x 4 lattice at the 
temperature range of 0.8 T/T pc ^ 1.4, where T pc is 
the pseudo-critical temperature on the lattice. 

Results of the lattice simulations have shown that mag- 
nitude of the gluon propagators increases at T < T pc . 
This can be regarded as contributions of the quark back 
reaction since the gluon propagators are almost insen- 
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TABLE I: Summary of fits with x 2 /Wdf 
types of fit ranges for each temperature. 


for several temperatures. We take two 


T/T pc 


p/Tpc | min 






d T z r4 x' 


! /JVdf 


0.82 


3.64 
4.64 


1.24(7) 0.52(1) 
1.60(8) 0.534(8) 


0.86 
1.13 


2.3(1) 0.827(9) 
2.1(1) 0.696(7) 


3.47 
1.41 


0.86 


3.82 
4.86 


2.1(2) 0.49(1) 
2.8(2) 0.46(1) 


1.73 
0.14 


1.82(5) 0.592(8) 
1.49(8) 0.460(8) 


3.51 
1.03 


0.94 


4.18 
5.32 


3.5(3) 0.47(1) 
4.3(3) 0.442(1) 


0.87 
0.13 


0.84(2) 0.287(8) 
0.75(1) 0.145(7) 


12.09 
2.90 


1.00 


4.44 
5.66 


4.8(4) 0.50(2) 
3.6(2) 0.33(1) 


1.70 
0.76 


0.61(1) 0.114(8) 
15.0(1.4) 0.59(1) 


11.28 
0.43 



TABLE II: Results of T' c evaluated 
from thermodynamic potential of the 
gluonic sector at several T/T pc . The 
first (second) parenthesis indicates 
statistical (systematic) errors. 



T/T pc 


T' c [GeV] 


0.82 


0.332(20)(10) 


0.86 


0.291(26)(11) 


0.94 


0.233(19)(24) 



sitive below the critical temperature in the case of the 
pure gauge theory [10]. This implies that there are many 
configurations which can couple to the gluon fields even 
at T < T pc since the gluon confinement is slightly bro- 
ken by the quark back reaction. The longitudinal (trans- 
verse) propagator has shown enhancement at T ~ T pc 
(1.2T pc ) and then both decrease when temperature in- 
creases. This means that not only the longitudinal propa- 
gator but also the transverse one should be of importance 
to describe the deconfinement phase transition since the 
enhancement of the propagators in momentum space in- 
dicates strong correlations in real space which may relate 
to divergence of the correlation length in the QCD phase 
transition. 

In order to estimate contributions of the quark back re- 
action to the deconfinement phase transition, we fit the 
propagators by "Gribov-Stingl" form and input these to 
gluonic potential obtained by the leading-order of the 
2PI formalism in the pure-gauge theory [14] . The critical 
temperature estimated in the gluonic potential decreases 
with temperature of the gluon propagators increasing at 
T Tp C This means that the enhancement of the quark 
back reaction reduces the critical temperature of the glu- 
onic sector. 

In this paper, we discuss qualitative contributions of 
the quark back reaction to the gluonic potential. It has 
however been found that the strong sensitivity to tem- 
perature even below T pc appears in two-flavor QCD, in 
contrast to the pure gauge theory. This implies that it 



is indispensable to introduce a temperature-dependent 
Gribov-Stingl form for a precise calculation of the criti- 
cal temperature in the gluonic sector T' c . In general, the 
temperature dependent parameters in the Gribov-Stingl 
form make the analytic calculation of the gluon potential 
(12) difficult. Then one needs to numerically calculate 
summations of the Matsubara frequency and subtract 
divergence of the propagators by using the same oper- 
ators at zero temperature on the lattice, which require 
huge numerical cost. Furthermore, to obtain abundant 
and detailed information of momentum dependence on 
the lattice, simulations on large lattice are important. In 
particular, it is known that the finite size effect of the 
gluon propagators becomes serious at infrared regions. 
Although we omit the data at small momentum, to draw 
a definite conclusion about the quark back reaction we 
need careful investigation with larger lattice size. We 
leave these to our future work. Establishment of reliable 
effective potential of the gluonic sector combined with the 
matter potential fi ma t will lead us to the understanding 
of the phase structure of QCD at finite temperature and 
density. 
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